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Abstract 

The significant discussion about the possible chaotic behavior of 
the mixmaster cosmological model due to Cornish and Levin [J.N. 
Cornish and J.J. Levin, Phys. Rev. Lett. 78 (1997) 998; Phys. Rev. 
D 55 (1997) 7489] is revisited. We improve their method by correcting 
nontrivial oversights that make their work inconclusive to precisely 
confirm their result: "The mixmaster universe is indeed chaotic". 
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The Bianchi IX (mixmaster) cosmological model was investigated by Be- 
linskii, Khalatnikov and Lifshitz (BKL), clarifying the nature of cosmological 
singularities (see and Refs. therein), and by Misner, in his attempt to solve 
the cosmological horizon problem . The subsequent emergence of chaos in 
maps associated with this model attracted the attention of several other re- 
searchers in the following years. Later, it was realized that standard chaotic 
indicators used in previous works, like Lyapunov exponents, are not invariant 
under space- time diffeomorphisms (see for a comprehensive discussion). 
Since then, a number of contributions looking for coordinate independent 
manifestations of chaos have appeared in the literature, and the mixmaster 
cosmology has become one of the main paradigms of deterministic chaos in 
General Relativity. Even though some important results have been obtained, 
a satisfactory comprehension of chaos in this system remains to be achieved 
(see 1^ for a critical review). Among the most promising proposals we cite 
the efforts of Szydlowski and collaborators in treating the problem in terms 
of (invariant) curvature P]. This approach, however, breaks down in the mix- 
master case 0]. Another significant progress was recently made by Cornish 
and Levin (CL) 0. Employing fractal techniques, the authors obtained evi- 
dences of chaotic transients in the mixmaster model. Their work is, however, 
inconclusive due to conceptual fiaws. Accordingly, the existence of invariant 
chaos in the mixmaster cosmological model is still an open question. 

In dynamical systems, chaos seems to be the rule rather than the excep- 
tion (at least for generic autonomous systems with three or more equations). 
Based on this fact, there is small risk in declaring in advance that a system is 
chaotic. In addition, even in nonchaotic systems, technical errors in general 
lead to some kind of numerical chaos. Therefore, the relevance of stating 
that a particular system is chaotic stays on the precision of the proof and 
the nature of the chaos revealed by it ||^. In the case of mixmaster universe 
it stays also on the need of a proper meaning for invariant chaos in General 
Relativity. 

The aim of this communication is to show unequivocally that the mix- 
master dynamics is chaotic in a precise sense. First we review critically the 
work of CL, pointing out some key oversights. Second, we propose modifi- 
cations and corrections of the CL's procedure. Next, this new procedure is 
employed to the approximated dynamics and then to the exact dynamics. 
From the full dynamics analysis we conclude that the mixmaster model does 
evolve chaotically. 
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The work of CL is based on enlarging the three outcomes of the mixmaster 
universe and in studying the invariant set through fractal methods of chaotic 
scattering. The existence of a fractal structure in the invariant set would 
imply chaotic transients. First they studied an approximation of the exact 
dynamics, the Farey map, whose invariant set could be exhibited explicitly. 
For this map, the topological entropy and the information dimension of the 
invariant set were computed. Then, by numerically computing the informa- 
tion dimension for a choice of the outcomes, the full dynamics was studied. 
Based on the noninteger value of the resulting information dimension, the 
authors declared that the system is chaotic. 

In checking the CL's papers we see that: (1) Concerning the Farey map, 
the invariant set exhibited, a Cartesian product of periodic points of the 
epoch-era BKL map is not invariant. In contrast to the authors' belief, 
only points with the same period in both Cartesian components are periodic 
points of the Farey map. Consequently, the conclusions obtained from this 
'invariant set' are wrong or at best suspicious. For instance, their result 
for the topological entropy is twice as bigger as the correct value because 
what they found was the square of the number of periodic points that really 
exist. (2) In the full dynamics counterpart, the information dimension was 
computed from the box-counting dimension of the basin boundaries, based 
on a result presented in Ref. - the information dimension of a set may be 
given by the box-counting dimension of an optimal f Taction of the set (see 0, 
p. 83) - that does not apply to this case because the optimal condition is not 
satisfied. Moreover, the initial conditions, given near the moment of maximal 
expansion, were parametrized by parameters from Chernoff-Barrow (CB) 
map and the outcomes were defined in terms of one of these parameters [p!0|] . 
Since the CB map cannot be applied far from the cosmological singularities, 
this parametrization is violated during the integration and is recovered only 
when the Kasner approximation becomes valid. Accordingly, the outcomes 
are not well defined all the time and they make sense only when a Kasner 
approximation is made. Consequently, the full dynamics analysis is not really 
full. (3) Both approximated and full analyses involve the idea of a full repellor 
as a set of the system itself. This concept is unclear because repellors are 
functions of the outcomes, and therefore are sets that are defined only after 
the definition of the outcomes. The origin of the problem seems to be in 
some confusion between nonattracting invariant sets (repellors) and the set 
of periodic points. 
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When considered the above oversights, the CL's work becomes inconclu- 
sive. In spite of that, the authors left a great contribution: The suggestion 
of using fractal invariant methods to investigate the mixmaster mo del0. In 
this communication we will follow this idea as close as possible of the CL's 
work in order to allow a direct comparison by the reader. But different from 
CL, we will choose the box-counting dimension as the invariant indicator of 
chaos 10. The topological entropy and the information dimension, employed 
by CL, are unsuitable indicators of chaos in the mixmaster system: The 
topological entropy is a function of the time parametrization because of its 
dependence on the period of the periodic points; the information dimension 
is difficult to be computed since it is difficult to verify the optimal condition 
mentioned above (we stress that CL were not able to obtain such dimension 
in the case of the full dynamics). The box-counting dimension, on the other 
hand, is diffeomorphism invariant and can be easily computed by applying 
the uncertainty exponent method 0. 

The Einstein equations expressed in terms of the three spatial scales of 
the model (gij) =diag(a^, 6^, c^) lead to the motion equations 

(Ina^)" = (62 _ c2)2 _ a\ et eye, (1) 

together with the constraint equation 

{\na^y{\nb^y + 2a%^ - + eye. = 0. (2) 

A prime denotes d/dr, where dt = abe dr and t is the cosmological time. 
It is convenient to introduce ai, 012 and 03 as a^"-' = In a, oi^^^ = ln6 and 
QjC'^) = Inc in increasing order of the time derivatives. That is, di < 02 < c^s, 
where a dot denotes djd^ and f2 is the time ln(a6c)/3. The dynamics near 



the cosmological singularity is approximated by the CB map |jTO[, obtained 

(3) 



from a parametrization of the surface of section = 0, > by 
ai = 3figi(M,t;), 02 = 3fig2(M, f ), as = 3i7g3(-u, t;). 



dl=3pi(M), d2 = 3p2(M), d3 = 3p3(M), 

where gi = l/(l+f +uf ), ^2 = ^Ki{v^uv) j {\^v^uv),q':>, = 6i^2{v+uv)/{l+v+ 
uv),pi = -u/{l+u + u^),p2 = 6^2{l + u)/{l+u+u'^) + 6^3{u+u'^)/{l+u+u'^) 



^Fractal techniques were also introduced by Demaret and De Rop jlll to study an 
approximation of Bianchi IX cosmology. 
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andp3 = 5k2(m + m^)/(1 + m + m^) +5k3(1 + m)/(1 + m + m^), for < < oo. 
The evolution of the parameters u and v is then determined by the Farey 
map, 



. _ \ {u — l,v + 1) if u > 1 (change of epoch) 

{u,v) - < _ (t, + i)-i) if M < 1 (change of era). ^ ' 



I — > 



The scale factors! (cti, a2, a^) = evolve as di < 0:2 < «3 

«i < '52 < ois, where (011,02, as) = {a^'^\ a^^\ a^''^) and k = 2 except in the 
transition to the last epoch of each era, where (ai, ce2, as) = {a^'^\ a^'^\ a^'^^) 
and K, becomes 3. The Farey map can be equivalently defined as, 

p, _s _ j (u — l,v + 1) if M > 2 (change of epoch) , , 

^[u,v) = <^ ((m - l)-\ v-^ + 1) if M < 2 (change of era), ^ ' 

ioT 1 < u,v < 00. Since F = f~^oFof {F and F are conjugate), where 
f{u,v) = [u — l,v — 1), the parametrization corresponding to F is not just 
(§) (as used by CL) but, (^ with qi and pi replaced by qi = qiof and 
pi =Piof, for i = 1,2,3. 

First let us consider the approximated dynamics in order to have some 
insight into the problem. In taking relative quantities in (^), the overall 
scale factor Q is canceled, and we are led to a stationary phase space u x 
V where three outcomes can be identified: (m, 03) = (00, a*^"-*), {u,as) = 
{00, a^''^) and (m, as) = (00,0^*^^). Since only rational values of u are led to 
these outcomes by forward iteration of the map, the corresponding invariant 
set contains almost every point of the phase space and it tells us nothing 
about chaos. However, some useful conclusions can be obtained if we enlarge 
the outcomes. It is made by defining the three outcomes for u > Uout — 1 
instead of m = 00, which transform the system in an exit (scattering) system 
(the original Farey exit map was defined by CL [^). These outcomes are 
equivalent to Q;3/(Q;i+d2 + tt3) > Pa (wout), where p3('u) = {u+u'^)/{l+u+u^). 
The dynamics of this exit system is determined by a nonattracting invariant 
set - the repellor [^j. The repellor is a zero measure set consisting of a 
countable number of unstable periodic orbits surrounded by an uncountable 
number of nonperiodic orbits. A fractional box-counting dimension for the 



^The indices /z, v, p are permutations of a, b, c. 
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repellor implies^ the exit system to be chaotic in the sense of having a chaotic 
transient: The system evolves chaotically for a period of time before being 
scattered. Concerning the original system, the system is typically never 
scattered (except for particular choices of the initial conditions) and hence it 
would evolve chaotically forever. Therefore, if the exit system is chaotic, the 
original system is also chaotic. 

In Fig. 1 we show, in three different colors, the basins of initial conditions 
for the outcomes defined by Uout = 7 (the same used by CL 0). It can 
be shown that magnifications of the basin boundary lead to figures with 
the same complex structure of Fig. 1, where the three colors (basins) are 
always present. The basin boundary is the future invariant set, and its 
box-counting dimension estimated from the uncertainty exponent method 
results Dq = 1.87±0.01. Since the repellor is the intersection of the future 
and the past invariant sets, it is equivalent to study any one of these sets. We 
prefer the future invariant set because of its easier geometric interpretation. 
Therefore, the above result for Do is enough to conclude that the original 
map is chaotic. In varying the outcomes, Dq becomes an increasing function 
of Uout as showed in Fig. 2a. When the outcomes are made arbitrarily 
small. Do approaches 2. This behavior is completely different in a nonchaotic 
system, since in that case Do is integer and does not approach the phase space 
dimension when the outcomes are reduced. 

Now we repeat the procedure for the full dynamics, numerically integrat- 
ing the mixmaster equations. In order to allow comparisons with the approx- 
imated study, we take the initial condition defined by (|^) for 1 < -u, < 2 
(k = 2), and Q the negative solution of the constraint equation. The integra- 
tion is performed in the negative direction of Q (approaching the singularity). 
Since u is not well defined during the integration, the best choice for the out- 
comes is d3/(di + a2 + as) > Psiuout)- In the limit Q — > — oo these outcomes 
correspond to those defined for the Farey map. The full dynamics is, how- 
ever, considered for a continuous time and the outcomes are defined not only 
on a surface of section. In Fig. 3 we show the basins of initial conditions for 
Uout = 7. Fig. 3 may be compared with Fig. 1. Small differences between 
these figures were expected, since the latter is an asymptotic approximation 

■^This implication is of typical validity and includes the extreme case, when the box- 
counting dimension of the repellor equals the phase space dimension. Pathological coun- 
terexamples can be constructed, but do not seem to represent realistic systems . 
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while the former is exact and refers to initial conditions taken near the maxi- 
mum of expansion. The complicated structure of Fig. 3 is also present in any 
magnification of the basin boundary - see Fig. 4, where we show a 100 times 
magnification of Fig. 3. The numerical computation of the box-counting di- 
mension corresponding to Uout = 7 results Dq = 1.87 ± 0.01, consistent with 
the approximated value. In fact, the differences between the approximated 
and full fractal dimensions are quite small, as can be seen in Fig. 2. More- 
over, the system has the so called Wada property [|13|: Every point on the 



boundary of a basin is also on the boundary of the other two basins. The 
corresponding fractal figures of CL do not present this property 0. 

From the above analysis follows that the mixmaster dynamics is in fact 
chaotic. This conclusion is independent of the particular space-time param- 
eters used because Dq is invariant under space-time diffeomorphisms. The 
numerical values obtained for Dq depend on the choices made for the out- 
comes. But, once found a set of outcomes for which the resulting Dq is 
fractional and goes to the phase space dimension in the limit "outcomes 

0" , the system is unequivocally chaotic. 

The authors thank R.R.D. Vilela and M.T. Ujevic for their valuable com- 
ments, and also Fapesp and CNPq for financial support. 



7 



References 

[1] V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, Adv. Phys. 19 (1970) 
525. 

[2] C.W. Misner, Phys. Rev. Lett. 22 (1969) 1071. 

[3] Deterministic Chaos in General Relativity, eds. D. Hobill, A. Burd and 
A. Coley (Plenum Press, New York, 1994). 

[4] A. Burd in Ref. [3]; G. Contopoulos, B. Grammaticos, and A. Ramani, 
J. Phys. A 28 (1995) 5313. 

[5] M. Szydlowski and A. Lapeta, Phys. Lett. A 148 (1990) 353; M. Szyd- 
lowski and A. Krawiec, Phys. Rev. D 53 (1996) 6893. 

[6] J.N. Cornish and J.J. Levin, Phys. Rev. Lett. 78 (1997) 998; Phys. Rev. 
D 55 (1997) 7489. 

[7] P. Holmes, Phys. Rep. 193 (1990) 137. 

[8] J.D. Farmer, E. Ott and J.A. Yorke, Physica D 7 (1983) 153. 

[9] E. Ott, Chaos in Dynamical Systems (Cambridge University Press, New- 
York, 1994). 

[10] D.F. Chernoff and J.D. Barrow, Phys. Rev. Lett. 50 (1983) 134; B.K. 
Berger, Phys. Rev. D 49 (1994) 1120. 

[11] J. Demaret and Y. De Rop, Phys. Lett. B 299 (1993) 223. 

[12] G. Troll, Chaos Soliton Fract 7 (1996) 1929. 

[13] D. Sweet, E. Ott and J.A. Yorke, Nature 399 (1999) 315. 



8 



Figure 1: Portrait of the basins of the approximated dynamics for Uout = 7. 
The initial conditions were chosen on a grid of 400 x 400, for ai = a^°'\ 
OL-i — and 0:3 = q;^'^\ Regions in gray, black and white correspond to 
orbits that escape for 6l2,I{pl\ + q;2 + ol-^) > 56/57 with as — a^°'\ — a^^^ 
and as — a^'^\ respectively. 



Figure 2: Box-counting dimension (Dq) of the basin boundary (future invari- 
ant set) as a function of the outcome parameter (uout) for: (a) Approximated 
dynamics (triangles); (b) Full dynamics (squares). The size of the triangles 
and squares corresponds approximately to the statistical uncertainty of Dq. 



Figure 3: Portrait of the basins of the full dynamics for Uout — 7. The initial 
conditions were chosen on a grid of 400 x 400, for «i = a^"'\ «2 = ct^'^^ and 
0:3 = a^^\ Regions in gray, black and white correspond to orbits that escape 
for a3/{ai + a2 + as) > 56/57 with as = a^"'\ as = a^''^ and as = a^'^\ 
respectively. 



Figure 4: A 100 times magnification of a portion of Fig. 3. 
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